Coherence Constraint on the Existence of Cosmic Defects 
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We devise a measure n{r) to quantify a robust feature in the coherent cosmic perturbations, 
namely the narrowness of the first peak in the power spectrum of the Cosmic Microwave Background 
(CMB) anisotropy. A maximum-likelihood analysis using the WMAP data shows that the power 
fraction of any defect models at the first peak is persistently less than 4.5% at the 68% confidence. 
Our approach and results are insensitive to the uncertainties in the theoretical predictions for defect 
models, and thus the most robust to date. We show that to convert these results into realistic 
constraints on theories such as infiation, string cosmology, and SUSY GUTs, a more robust study 
for the cosmic-string-induced CMB around the degree scale is required. 

PACS numbers: 98.70.Vc, 98.80.Cq, 98.80.Es 



A. Introduction: Cosmic defects have been an in- 
teresting subject for different reasons through different 
decades. They were first introduced around the 70's as 
an inevitable consequence of various theories such as the 
quantum field theory and the unification theory P, 01 • 
Then motivated by the anisotropy in the Cosmic Mi- 
crowave Background (CMB) observed by COBE in 1992 
, they competed with inflation as the dominant mech- 
anism for the formation of structures in the universe 
H H 0, 0,11 0. They started losing this contest in 
the recent years when the flood of high-precision data 
showed great agreement with inflation (e.g., 0, lll|). 
Nevertheless, defects have been still playing a key role in 
testing some of the fundamental physics, especially those 
related to inflation and string cosmology. In general, for 
example, various defects may form as hybrid inflation 
ends • In string theory, colliding branes can not only 
bring the brane inflation to the end but also potentially 
generate cosmic strings [isl . In the Supersymmetric 
(SUSY) Grand Unified Theories (GUTs), which embed 
the standard model into the string theory, cosmic strings 
can also be formed at the end of hybrid inflation !l5l. llq . 
Therefore, we can use observational data to place limits 
on the energy density of cosmic defects, so as to confine 
the energy scales of the associated physical mechanisms. 
In this paper, we will employ the CMB data for this pur- 
pose. 

There are two main types of defects: global defects 
(Models I-IV in Tab.P) and local cosmic strings (Models 
V-VII). Their resulting CMB power spectra are shown 
in Fig. n] The CMB power spectrum is defined as 
Ce = (|a^mp), where a^m's are the multipole-expansion 
coefficients of the CMB temperature anisotropy, and £ is 
the multipole number. By comparing these results side 
by side, we find that the results for the 0(4) textures 
(Model III) in Ref.0 and 's'l agree weU at the 10% level 
for £ < 1000. On the other hand, the uncertainty in the 
theoretical predictions for cosmic strings remains large. 
Models VhJ and VM are both based on a toy model 



using filaments to mimic strings, while Model VIl'^ uses 
string simulations extrapolated from a non-expanding 
universe. These predictions disagree not only in shape 
but also in amplitude. For example, although based on 
the same approach, the predicted Ci of Model VI 6] is 
higher than that of Model V 0| by a factor of 4 and 16 
at £ = 10 and 700 respectively (for the same string linear 
energy density fi, where Ce oc ii"^; see Fig.QJ. Depending 
on the velocity of the decay products of the long strings, 
the amplitude and the shape of Ce may vary significantly 
(see Model VII in Fig. nj. These theoretical uncertain- 
ties have been overlooked in most of the recent studies 
about the observational constraint on and other related 
model parameters (e.g., ll3,0,^3)i dramatically weak- 
ening the conclusions therein (e.g., < 2 x 10~^ at 
99% confidence in Ref.[ilEl, which used WMAP and 
SDSS data). Note that work employing full string sim- 
ulations in an expanding universe has also been pursued 
d 0, m US 113 but onlyfor £ up to a few tens 

due to the numerical limitation |22l |. This £ range is not 
suitable for us to investigate the observational first peak 
at « 220 (see later) , so we will not consider this model 
here. It is still worthy mentioning that the Cio in Ref.|2^ 
is about 8 times higher than that in Ref. for the same 

In face of the above uncertainties in the theoretical 
work, we devise a new strategy to constrain the en- 
ergy level and thus the existence of defects. This strat- 
egy focuses on a robust feature of the coherent cosmic 
perturbations — the narrowness of the first peak in Ce. 
B. Quantifying the coherence: We first define Ce — 
£{£+l)Ce/2Tr. If the apex of the first peak in Ce is located 
at £p with the height Cp = C{£p), and the width of the 
peak at the height rCp (0 < r < 1) is A£{r) (see Fig.EJ, 
then the new measure is defined as 



n{r) 



A£{r) 



(1) 



As will be justified and explained later, the shape of the 
function n[r) is determined by the nature of whether the 
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FIG. 1: CMB power spectra of global defects (blue dashed; 
Models I-IV downwards), local cosmic strings (red dot- 
dashed, red dotted, and red thin solid), the concordance 
model (black thick), and the WMAP (boxes; at l-cr). The 
19 empty boxes are used for the analysis of n(r) and /22o- 
Model V is normalized with = 2 x 10"*'; Models VI and 
VII are both with Gfi = 7 x 10~^ for the ease of comparison. 



TABLE I: Cosmic defect models investigated here, and their 
power fraction /f2o a.t £ — 220, their temperature anisotropy 
at ^ = 220 and 10 (in ^K), as constrained by WMAP, 
all at the 68% confidence. 



X 


Models 


/220 


V 220 


1 fx 

V 10 


constraints 


I 


strings [5] 


< 4.5% 


< 16 


< 13 




II 


monopoles[^ 


< 4.4% 


< 16 


< 14 




III 


textures [^,^ 


< 3.7% 


< 14 


< 15 


e < 8.2 X 10"'' 


IV 


0(6) text.ia 


< 3.0% 


< 13 


< 14 




V 


loc. COS. str.fl7] 


< 3.4% 


< 14 


< 10 


Gfi<7 X 10"'^ 


VI 


loc. COS. str.f6] 


< 2.1% 


< 11 


< 4.6 


Gfi<2x 10^'^ 


VII 


loc. COS. str.fQ] 











CMB perturbations are coherent or not. It depends only 
weakly on the geometry of the universe, because a change 
in geometry will stretch both A£ and £p in the same way, 
so that their ratio n{r) remains very much unchanged. 

To justify the robustness of this coherence feature in 
n(r), we consider the inflationary models for a reasonably 
large region in the parameter space: < f^A ^ 1 (energy 
density parameter for cosmological constant) , 0.1 < < 
1 (for cold dark matter), 0.01 < Jib < 0.2 (for baryons), 
0.5 < /i < 0.9 {Ho = lOOft. km s~iMpc"i where Hq is 
the Hubble parameter today), 0.8 < rig < 1.2 (scalar 
spectral index), < t < 0.25 (optical depth). Note 
that this gives 0.11 < < 2.2 (total energy density 
parameter). First we compute all the CMB power spectra 
|24| within this parameter range, and then the resulting 
n(r) for each spectrum. These n(r)'s occupy the blue 
shaded region in Fig. |21 It is clear that the coherence 
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FIG. 2: Illustration for the definition of n(r) in Eq. Q. 
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FIG. 3: The n(r) of the concordance model and the defects 
(same line styles as in Fig. . Additionally, the inflationary 
models spam over the blue shaded region, and the red lined 
regions are the 68% and 95% C.R.'s from WMAP. 



nature of the inflationary models confines the n{r) within 
a very limited region on the r-n plane. 

On the other hand, we compute the n(r) for the defect 
models in Tab. Obviously the peaks in Ci for these 
incoherent models are much broader (Fig. resulting 
in an n{r) of at least twice higher at a given r (Fig.|3Jl. 
A comparison between these lines and the shaded blue 
region shows the capability of n(r) in discriminating the 
coherent and the incoherent models. This discrimination 
is independent of the cosmological parameters. 

Theoretically this coherence feature in the n{r) can be 
manifested by the following intuitive approach. In the 
synchronous gauge for a perturbation theory lisf . 
it is straight forward to show that with inflation the 
power spectrum of photon energy contrasts at the last- 
scattering epoch 771s can be approximated by the form: 



P{k) oc A:"= 



cos 



1 + [kmsY 



(2) 



where k is the comoving wave number. The exponential 
term accounts for the photon-diffusion damping, with rj^ 
typically one order below -qis jl^ . The cosine oscillations 
are due to the coherence nature of the perturbations, 
resulting from the fact that the inflationary perturba- 
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tions on the same scale {k~^) entered the horizon at the 
same time k~^) and thus oscillated coherently. On 
the other hand, for models with only defects, the per- 
turbations on the same scale were seeded over a range 
of time after the horizon crossing, and thus oscillated in- 
coherently. These incoherent oscillations effectively re- 
move the cosine dependence in Eq. leaving only 
a broad peak around the scale rjis- We then compute 
Ci (X J P{k)ji{kr]o)k'^dk, where ijq is the comoving ra- 
dius of the last-scattering sphere and ji is the spherical 
Bessel function. We find that the resulting n{r) for the 
coherent case is slightly lower than the blue shaded region 
in Fig. 131 due to the omission of the Doppler effect and 
the Integrated Sachs- Wolfe effect. These two effects are 
sub-dominant in the first peak and thus in n(r). It is im- 
portant to note that the cosine dependence in P{k) gives 
not only the periodic peaks in Ci as commonly known, 
but also the narrowness of the first peak as we try to 
advocate here. On the other hand, the incoherent case 
leads to a much broader peak in CV, resulting in an n{r) 
of at least twice larger than the coherent case. 

The cosmic perturbations may be attributed to both 
coherent and incoherent mechanisms. Given the fact that 
the former is the dominant, an increasing fraction of the 
latter in a mix model will broaden the width of the first 
peak, resulting in a higher value of n(r). In the following, 
we will invoke this feature to constrain the fraction of the 
incoherent perturbations and thus the defects. 
C. Observational constraint: From observations, we 
can estimate n{r). Among the large amount of CMB 
results that measure the first peak in Cg, we choose the 
WMAP because it is the most stringent to date. We 
will denote its result as (^^(i), €i(i.),Ci, tJi), where i indi- 
cates the i-th £ bin, and denote respectively the 
smallest and the largest £ within the bin, Ci is the binned 
power, and cr^ is the 1-a error in Ci. For a given point 
A{£'p,Cp) on the £-Cg plane, the probability that it is the 
apex of the first peak is 

Pi(A) =pd/(C;;Ca,a,)[]cd/(C;;C„aO, (3) 

where a denotes the a-th £ bin within which £p lies (i.e. 

G [^a(i)) C(r)])) pdf(x; /i, a) is a normal probability dis- 
tribution function with mean /i and standard deviation 
(7, cdf{x;fi,a) is a normal cumulative probability distri- 
bution function. Then at the height rCp (0 < r < 1), the 
probablility for point B{£[,rCp) being the left-end point 
of the peak and A£' being the width of the peak is 



P2iB,A£'\A)^ II cdf{rC;;C,,a, 

i,(,)<e[ or £;+A^'<£i(i) 

n [l~cdf{rC'p;C,,a,)] 



X 



where b and c denote respectively the £ bins within 
which the £[ and the £[ + A£' are, and /(a:; rCp, •^i(i), ^^(2)) 
is a linear-interpolation function between — 
cdf{rCp;Ci,(Ti)) and {£i(2), cdf{rCp]Ci, ai)) according to 
X. Here the condition £[ < £p < £[ + A£' is also imposed. 
Finally, given the observational data, the probability for 
a point on the r-n plane to be present is 



P{r, n) 



E 



Pi{A)P2{B,A^'\A). 



(5) 



For each r, we compute the 68% and the 95% confidence 
regions (C.R.) for n, as shown in Fig. 01 Note that for 
the WMAP data, we use only the 19 bins (i = 12-30; 

< £ < 379) that best confine the first peak. These 
bins are shown as the empty boxes in Fig. ^ 

It is quite clear in Fig. O that the observational data 
strongly favor the inflationary models, regardless of the 
geometry of the universe. Thus a defect-dominated uni- 
verse with a very closed geometry, which serves to shift 
the first peak from several hundreds to around 220, is to- 
tally inconsistent with the observations. Also plotted in 
Fig. 131 is the n{r) of the concordance model based on the 
data of WMAP, CBI, ACBAR, 2dF, Ly-alpha Its 
power spectrum is shown in Fig. ^ Its n(r) is roughly 
within the 68% C.R. of the data. We also note from 
Fig. 131 that on average the data favor a slightly narrower 
peak than what is predicted by the inflationary models. 

To see how stringent the observed narrowness of the 
first peak can constrain the existence of cosmic defects, 
we further consider the following mixed power spectrum: 



nt 



A 



^220 



I-VI, (6) 



xl{£[; rC;, 4(i), 4(r))/(^; + A^'; rC'^, 4(r), 4(i)), (4) 



where the superscripts 'inf and 'x' denote the contribu- 
tion from inflation and from the defect models in Tab. jH 
respectively. Model VII is not considered here because 
Ref.jil does not provide information about the depen- 
dence of the Ce on the cosmological parameters. Here 
/220 is nothing but the power fraction of the defects at 
the first peak {£p « 220), and A is an arbitrary normal- 
ization factor. We use the WMAP data to constrain /220J 
again using only the 19 £ bins. A standard maximum- 
likelihood analysis is performed for the parameters of 
/22O' ^A, f^c, f^b, h, Us, r, with the flat-geometry 
constraint that ilt = Oa -I- fie + f^b = 1. We then 
marginalize over all the parameters except for /f20j leav- 
ing only a one-dimensional likelihood C-{f 220)- At the 
68% confidence level (C.L.), the estimated values of /f2o 
are shown in Tab. Ql Two robust features in all results 
for different models are that the likelihood C{f22o) peaks 
below /I20 = 1.2%, and that /220 < 4.5% (or equiva- 
lently ^Jc^ < 16^K) at the 68% C.L. This means that 
the defect models contribute at most 4.5%, preferably be- 
low 1.2%, in the observed C220: or equivalently that their 
contribution to the CMB anisotropy around the degree 
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scale is no more than 16/iK, preferably below 8.2/iK. The 
y^C^Q in Tab. n is extrapolated from the value of \/C^ 
taking the corresponding model shape ofCf. We see that 
■y/Cfp < 15/iK (around the COBE scale) for all models. 
D. Discussion and Conclusion: We emphasize that 
although the shape and the amplitude of vary a lot 
among different defect models, the observational con- 
straint on their /220 remains persistently below 4.5%. 
When we replace the in Eq. ijHJl with a Ce which has 
dCe/d£ = m, the same conclusion persists for — 10 < 
m < 10. This indicates that any further improvement 
on the detailed predictions of Cf for defect models can 
only marginally affect our results here. In theory this 
is naturally due to the fact that the predicted Cg^s are 
lack of curvature around the observational first peak at 
£ « 220, so that a larger /fjo will result in a higher n{r), 
which potentially violates the observational data. 

Our results can be converted to place limits on the 
parameters in various theories. For example, the re- 
sult 



< 15/J,K can be converted to e = 47rG0Q < 
8 X 10^^, leading to a constraint on the energy level of 
symmetry breaking 0o < 9.9 x lO^^GeV. This is consis- 
tent with, though more conservative than, the result in 
Ref.ji^, where the full WMAP data were used to yield 
/{" < 13%, or equivalently y/C^ < lO/xK, at 95% C.L. 

For local cosmic strings, naively we can take the more 
conservative result for fi, from Model V (see Tab. P), 
to yield a constraint on the symmetry-breaking scale 
77 ~ fj}^^ < lO^^GeV. In the context of string theory, 
where the strings are D-branes, it places an upper limit 
onto the superstring scale ^ (/x/2)i/2 < 7 x lO^^GeV. 
In the context of SUSY GUTs, a combination of our fur- 
ther deduced result Q^ms-ps = (^^)q ^ 4.4/xK and 
the study in Ref.[l6l| yields an upper limit on the mass 
scale of the F-term inflation A/ < 2 x lO^^GeV, and 
on the superpotential coupling k < 2 x 10^^. For the 
D-term inflation, we get constraint on the gauge cou- 
pling 5 < 1.4 X 10~^, and on the superpotential couphng 
A < 2.2 X 10^^. We emphasize that all these constraints 
still carry considerable uncertainties due to the large un- 
certainties in the predicted C220 for strings, as addressed 
in the introduction. This uncertainty is propagated into 
the constraint on /.t and thus on the deduced parame- 
ters. This is readily seen by the inconsistency between 
G/i < 7 X 10-'^ and 2 x 10"'^ for Models V and VI respec- 
tively (see Tab.nj. Therefore a more robust predictions 
for strings around the degree scale will dramatically set- 
tle this uncertainty. For the same reason, whether or not 
the values of G/i = 5.9 x 10"'^ Hi] and 4 x lO"'^ p3| that 
were required to explain the recent anomalous observa- 
tions in the gravitational lens systems is consistent with 
the current CMB constraint remains an open question. 



models at the first peak is persistently less than 4.5% at 
the 68% C.L., in spite of the large uncertainties in the 
current theoretical study (see Tab.lJfor the more detailed 
results). Hence our results should be the most objective 
and thus the most robust to date. 

We acknowledge the support from the National Science 
Council of Taiwan (NSC 92-2112-M-002-047; NSC 93- 
2112-M-002-015). 



Nevertheless, we are confident that a/CI^ — 14/iK. 

In conclusion, we invoke a coherence feature — the nar- 
rowness of the first peak in Ci — to constrain the existence 
of defects. We show that the power fraction of any defect 
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